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The one-loop dynamical polarization function of graphene in an external magnetic field is cal- 
culated as a function of wavevector and frequency at finite chemical potential, temperature, band 
gap, and width of Landau levels. The exact analytic result is given in terms of digamma functions 
and generalized Laguerre polynomials, and has the form of double sum over Landau levels. Various 
limits (static, clean, etc) are discussed. The Thomas-Fermi inverse length qF of screening of the 
Coulomb potential is found to be an oscillating function of a magnetic field and a chemical potential. 
At zero temperature and scattering rate, it vanishes when the Fermi level lies between the Landau 
levels. 

I. INTRODUCTION 

The fabrication of graphene [l| initiated extensive theoretical and experimental studies of its remarkable electronic 
properties aimed at promised applications of this material in next-generation electronic devices. The non-interacting 
charge carriers in single layer graphene are described by the analogue of the Dirac equation for the massless fermions 
with the relativistic-like linear spectrum Q and a vanishing density of states at zero doping. In the presence of 
the external magnetic field the spectrum of these Dirac quasiparticles has the form of relativistic Landau levels, in 
contrast to the equidistantly spaced levels in a usual two-dimensional electron gas. These peculiar features of the non- 
interacting charge carriers in graphene result in several interesting physical phenomena such as the unconventional 
quantum Hall effect , the universal optical conductivity @, Q and magneto-spectroscopy || d, ■ 

Although these and other electronic and transport phenomena in graphene are well described in terms of free Dirac 
quasiparticles, the effects of interactions, in particular, the Coulomb interaction, are not settled yet. The vanishing 
density of states at the Dirac point ensures that the Coulomb interaction between the electrons remains unscreened 
due to vanishing of the static polarization for q — » The large value of the unscreened coupling constant 

g = e 2 /Hvf, where e is the electron charge, vf ~ 10 6 m/s is the Fermi velocity, could lead to instability in pristine 
graphene and formation of excitonic condensate and a quasiparticle gap, followed by quantum phase transition to an 
insulating phase above some critical g c . This possibility is studied in a series of theoretical works [lU, [p| (see, also 
recent papers [lij]) but experimental evidence for such an insulating phase is still absent [l5| . 

The screening of Coulomb potential due to the many-body interactions is determined by the polarization function 
which is also an important physical quantity for the spectrum of collective excitations (plasmons). This function in 
monolayer graphene without a magnetic field has been studied in one- loop approximation in Refs. (T3, IT5, fl7j ■ In the 
presence of an external magnetic field, it was calculated in [l8[ at zero temperature and impurity rate with the result 
given by the double sum over the Landau levels. The similar expression was also obtained later in (l9j . where it was 
employed to study the spectrum of collective excitations in a magnetic field. However, to the best of our knowledge, 
the most general expression for the dynamical polarization in the presence of finite temperature, chemical potential, 
impurity rate, quasiparticle gap and a magnetic field was not given in the literature. 

The present paper deals with this more general case. The paper is organized as follows. In Sec. [II] we describe the 
model used and present our main result for the polarization function. We consider the clean graphene limit of this 
function in Sec. IIIH In Sec. [TV] we focus on the static screening properties of graphene. Then, in Sec. El we discuss 
some other limits of the polarization function, and in Sec. IVII we give the brief summary of our results. Finally, we 
provide the details of the calculations in the appendices A and B. In the appendix A we derive the expression for the 
dynamical polarization as a double sum over the Landau levels while in the appendix B we employ the Schwingcr 
proper time method to get a double integral representation for the polarization. 
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II. MODEL AND GENERAL EXPRESSION FOR POLARIZATION FUNCTION 



The Lagrangian describing the non-interacting Dirac quasiparticles confined to the graphene plane, in an external 
magnetic field, reads (we use the units h = c = 1) 



V(^t - i/x) + ivFl{V + ieA cxt ) - A , (1) 



=i 



where = {^^a^kb^k'b^k 1 a) ^ s ^ ne four-component wave function describing the Bloch states on the A and 
B sublattices and in the vicinity of K and K' points in the momentum space. = ^^.7° is the Dirac conjugated 
spinor, a is the spin variable, and gamma-matrices 7^ = 03 (8 (03, 202, form the reducible 4x4 representation 

in 2 + 1 dimensions. 

We will neglect the Zccman splitting which in graphene is very small (~ 1.34i?[T] K) compared to the distance 
between the zeroth and the first Landau levels (~ 424-^/B[T] K). Therefore, the electron spin results in only the 
degeneracy factor (number of flavors) Nf = 2. We have also included the gap term A which can be induced in 
graphene by placing it on a top of an appropriate substrate [2p| that breaks the sublattice symmetry, or can be 
generated dynamically in magnetic field (the phenomenon of magnetic catalysis) [l2l . [l3j . The external magnetic field 
B = V x A oxt is applied normally to the graphene plane and the vector potential is taken in the symmetric gauge 
A ext = (—By/2,Bx/2). The chemical potential [i can be varied by applying the gate voltage. 

The Green's function of Dirac quasiparticles described by this Lagrangian in an external magnetic field reads 

G(t - t', r; r') = cxp(-ierA cxt (r'))S{t -t',r- r') , (2) 



where S(t—t', r— r') is the translation invariant part of the propagator. Using the expression for S(icu s , q) from [13ll2l| 
we obtain for the propagator in the configuration space (in the Matsubara representation) 



qr , ( r 2 x ^ [ 7 °(io; m + /1 + iT n sgn^ m ) + A]/f(r) + / 2 "(r) 

5^ m ,r) = — exp(--j^ (^ m + /i + i r n sgn Wm ) 2 -A^ ' "m = fa» + W , (3) 



n=0 



where T is the temperature (we use fee = 1), M n = yj2nv 2 F jl 2 + A 2 , and E n = ±M n are the energies of the relativistic 
Landau levels, 1=1/ y/\eB\ is the magnetic length. The functions /™2( r ) are defined as 

A"(r) = P.L n (^) + P + L n ^ (^) , P± = ±(l± i 7 V sgn(eB)) (4) 

= -y^TjlU^). (5) 

where L%(z) are the generalized Laguerre polynomials (by definition, L„(z) = L^(z) and L" 1 (z) = 0). 

The finite parameter r n has the meaning of the width of Landau levels or, equivalently, the scattering rate of 
Dirac quasiparticles. It is expressed through the retarded fermion self energy and in general depends on the energy, 
temperature, magnetic field, and the Landau level index. In our calculations, we are assume that the width is 
independent of the energy (frequency). 

The dynamical polarization determines many physically interesting properties, such as the effective electron-electron 
interaction, the Friedel oscillations and the spectrum of collective modes. The retarded one-loop dynamical polariza- 
tion function is given by the expression 

00 

U(in s ,q)^e 2 TN f d 2 re-^ r £ tr [j°S(iuj m , r) 7 °S(iu; m - -r)] , O s =27rsT, (6) 

m— — 00 

analytically continued from Matsubara frequencies to real f2 axis. Note that our definition of the polarization function 
differs in the factor of — e 2 from that used in Rcfs.[l8|, [l9j]. Details of the calculation of this function are given in the 
appendix, here we reproduce only the final expression and then analyze various limiting cases. Thus our main result 
reads 



2 AT 00 

n,n'— A, A' — ± 



z%',(n,r,»,T) 



\M n - \'M n , -a- i(r„ - r„o 
z-£'- A (fi,r,-/i,T) z^(n,r,n,T) + z-, x n r x (n,r,-»,Ty 



\M n - \'M„, - n- i(v n , - r„) XM n - \'M n , - n - i(r n + r„,) 



3 



where we have introduced the following notations: 



1 



and 



?^:(j/,A) = e-V"-"' 1 ' 



1 /i - XM n + + iT n 
' 1 2 2i?rT 



'1 /x - A'M n > + iT n , 
' 1 2 2i7rT 



(8) 



AY A 2 \ / n<! 
M n M n .) \ n> \ 

r 

|n-n L 



4T"''(y)l +(i-*on<) 



(n< ~ 1)! 
(n> - 1)! 



(9) 



+ /2M„M„- (n> -l)!^'- 1 (2/)I/ "< " (y) 



where y = I q /2, n> = max(n, n'), n< = min(n,n'), and V'C^) is the digamma function. Note the symmetry proper- 



ties of the function Q%n> (y> A) with respect to the exchange of indices A, A' and n, n' and Q^n' {Vi A) = Qnn' X (v> A) 
For the gapless graphene (with A = 0) the function © reduces to 



— y \n—n 

e <V 



\l + XX'6 0n> ) n< \ ln _ 



{n< - 1)! T \n-n'\, 



!T<~ n '(y) + aa'(i - <w<)W ^r^4"<"-i (f) 



Taking the limit of zero temperature, the expression (jHJ simplifies to 



Z£,(fi,I>,0) = — In 



/i - XM n + fi + iT r 
fi - X'M n ' + iT n i 



(10) 



(11) 



The polarization function ([7]) is an analytic function of fi without singularities in the whole upper complex half-plane. 
It depends only on the absolute value of chemical potential (that can be verified by the replacement A <H> —A', n <-> n') 
and obeys the relation II(— 57, q) = [II (fi, q)]* (can be verified by the replacement A -H- A' , n <H> n' and taking into 

account [Z^,(Q,,T, fi, t)]* = —Z~ n ,'~ (— fi, T, —^i,t)). At finite scattering rate, the polarization function ([7]) receives 
the contributions both from the inter- (with An ^ X'n') and the intra-Landau level (An = X'n') transitions. Note 
that Qqq~ X (u, A) = which reflects the fact that the levels with energies ±A belong to the different valleys, and the 
intcrvalley transitions are not incorporated in our model. 



III. CLEAN GRAPHENE 



In the absence of scattering of Dirac quasiparticles (T n = 0) the general expression (JJJ for the polarization function 
reduces by means of Eq. (|A17[) to the following form: 



n(fi,q) 



4tt/ 2 



E E 3w(2/,A) 



n,n'— A, A'— ± 



n F {XM n )-n F (X'M n ,) 
XM n - X'M n ' + fi + i0 



(12) 



where np(x) = [e^ x ~^/ T + 1] — 1 . One can easily see from the above expression that only the terms with An- ^ X'n' 
(corresponding to the inter-Landau level transitions) survive in the clean limit. However, this is not the case when 
the limit T — »• is taken after setting fi = (see Eq. ([T7| below). When both scattering rate and temperature are 
zero, it simplifies further to (the order of taking limits r„ — > and T — > is not important) 



n(fi,q) = 



e 2 N 



f 



E E 



)-+ 



;(y,A) 



e 2 i\r 

2 ^ n M n + M n , + C(fi + iO) + 4^T 2 

n,n'=0 C — i 



oo Np- 



f^ 2 -A 2 )EE E 



n=0n'=0A,C=± 

where we used the symmetry of the function Q^ n i{y, A) with respect to upper indices, 

"(/i 2 - A 2 )/ 2 



XM n - M n , + C(fi + iO) 



, (13) 



2v% 



(14) 



is the number of the highest filled Landau level (square brackets here denote the integer part of expression). For // < 
it is a positive number meaning the highest empty Landau level in the valence band. The first term in Eq. (|13l) describes 
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vacuum contribution and takes into account only interband processes while the second one represents intraband and 
intcrband contributions when the chemical potential lies in the conduction or valence band. Notice that this second 
term does not receive contribution from the terms with n = n' . A = +1. 
In the gapless case (A = 0) we have 

V ' 4; 4tt/ 2 ^ ^ M n + M n . + C(n + iO) 4ttZ 2 Z^ Z^ Z^ \M n - M n , + C(Q + i0) ' V ' 

n,n<=0(=± ^ y ' n=0n' = lA,C=± 

Expressions (flU|) . (|15p . coincide with the polarization function calculated in [3]. Refs. [l9| considered only gapless 
case and obtained the expression similar to Eq. ([15[) but with twice larger contribution of the lowest Landau level 
(n = 0), while the results of the papers [H, [23[ are completely different from ours. 

The static clean limit of the polarization function essentially depends on the order of taking limits 17 — > and 
r„ — > 0. Indeed, taking first the limit Q = 0, the expression for the polarization function (JT]) reduces to 

e2N f v^v-^aa/ a \ d ,,f l , M-AM n + «r. 



n=0 A=± v 

+ 4^2 Z. Z, «„n'ltf.^Jm AM„ — X'M n ' — i(T n — T n t) 

n,n —0 A, A — ± 
An ^A'n' 

where we took into account that the numerator of the third term in square brackets in Eq.([7]) vanishes at ft = 0. Here 
we also introduced the ultraviolet cutoff n c due to the divergence of the sum over the Landau levels at finite width 
r„. This cutoff is estimated to be n c ~ 10 4 /-B[T] due to finiteness of the bandwidth [l9| . The expression (JTS]) for 
static polarization is obviously a real function. In the clean graphene limit, L„ = 0, we get 

n( °' q) - 16^T 2-j Z, cosh 2(^y - 4^" Z. 2. Ow(».A) _ (17) 

n=0 A=± V 2T / n,n =0 A,A'=± 

An ^A n 

(the sum over the Landau levels is convergent). On the other hand, if we take limit f2 — > in (|T2|) . i.e., after 
setting r„ = 0, we obtain the expression (|17[) without the first term. This term gives the contribution from the 
intra-level transitions (n <-> n.) even at zero width of Landau levels. At zero temperature it turns into the sequence of 
delta-functions <5(/i ± M n ) and does not contribute at integer filling factors of Landau levels (y = 0, ±2, ±6, ±10, . . .). 
Therefore, for T — 0, we arrive at the same expression (Q~3|) with = 0. 



IV. STATIC SCREENING 

The screening of the static Coulomb potential <pQ (r) = Ze jr is determined by the static polarization function, 

oo 

Ze f d 2 q exp(iqr) Ze f dqqj$(qr) 



e Q J 2ir q + (27r/e )n(0, q) e J q + (27r/e )II(0, q) ' 



(18) 



o 



where Jq(z) is the Bessel function and £o is the background dielectric constant due to the substrate. 

In what follows we assume that even in the case of clean graphene, the limit L n — > of the polarization function 
is taken after the limit ft — > when calculating the screened potential (fT5|) . This order of limits, which leads to the 
expression (fTT)) for 11(0, q), seems to be more natural due to the fact that real graphene samples can not be completely 
free from impurities and some broadening of the Landau levels always occurs. 

In general case, the static polarization function is given by the expression (|16[) which does not have singularities 
(like, for example, the discontinuity of the second derivative at q = 2fj,/vp in the absence of magnetic field [HI). 
Therefore, the asymptotical behavior of the screened potential at small or large distances is determined solely by the 
asymptotics of 11(0, q) at large or small wavevectors, respectively. At large momenta we have the zero magnetic field 
result, 

e 2 iV f lal 

H(0,q)^^Z^, q->oo, (19) 



5 



and (|18j) implies 

<Kr)~— , r^O, (20) 

where 

ire 2 N{ 

e* = £ + — — *-« eo + 3.4 (21) 

is the "effective" background dielectric constant (iV/ = 2). At small values of a wavevector (q — > 0), the static 
polarization function (p~6]> behaves as 

11(0, q)~ ^-(q F + aq 2 ) , q^O. (22) 
In the case q F ^ 0, we find from (Ti"g)) the following asymptotical behavior 

0(r) ~ — ^ , r -> oo , (23) 

that describes the Thomas-Fermi screening in graphene (l6| . In contrast to the three-dimensional case where for 
nonzero charge density the Coulomb potential 1/r is replaced by an exponential decreasing one, in two-dimensional 
case we have 1/r 3 behavior at large r, which is the well known fact [24|. The strength of the screening is determined 
by the magnitude of q F = (27r/eo)n(0, 0). 

At zero momentum q = only the first term in Eq. (|16|) contributes and we get 

n=0A=± V 7 

The above polarization function obeys the following relation (2f| 

oo 

n(0,0) = e 2 A p( ^ T) = e2 / (25) 
OfJ. J 4Tcoslr (%£) 

— oo 

where D(e) is the density of states in graphene with impurities in magnetic field [26j |. 

D ( e ) = N f VV (2-S 0n )T n 

and p(/x,T) is the density of Dirac quasiparticles, 

oo 

p(M,T)= y del>(e) [«*■(«:) -0(-e)] . (27) 

— oo 

At zero temperature and finite scattering rate the quantity 11(0, 0) is proportional to the density of states at the Fermi 
surface, 



It is an oscillating function of chemical potential and a magnetic field [26|, and therefore, the screened potential at 
large distances oscillates with changing (i at a fixed magnetic field, or with changing B at fixed [i. 
For r„ — and finite temperature, (|24|) reduces to the expression 



km>- -SI EE^^e,. <») 

n=0 A=± V 2T / 
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-6 -4 -2 2 4 6 
FIG. 1: Long wavelength limit of the static polarization function 11(0,0) at Y — 0, A = 0, T — 0.08vf/1. 
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FIG. 2: Fourier transform of the screened Coulomb potential at F = 0, A = 0, T — 0-Olwjr/Z. Dot-dashed (black) line: 
(j, — 0.01vf/1, solid (red) line: /i = 0.1vf/1- Dashed (blue) line shows the unscreened case. 



which has qualitatively the similar oscillatory behavior, see Fig. Q] The weak magnetic field limit (I — > oo) of the 
above expression can be obtained by replacing n — > k 2 l 2 /2, with the sum turning into the integral over k, resulting in 



11(0,0) 



e 2 N f T 



In 2 cosh 



A 



2T 



A /IA + fl\ . 

—— tanh — — — + (u — > —Lb) 



•A 



(30) 



which agrees with [13l ]. 

Some numerical results for the screened Coulomb potential in the case of clean gapless graphenc at finite temperature 
are shown in Figs. [5J [3] (we used Eq = 1). The figure [21 shows the Fourier transform of the potential (pi 



4>(q) = 



2irZe 



S q + 2irn(0,q) ' 



(31) 



and the figure [3] represents the screened potential (fT5)) itself. While the asymptotics of 4>{r) is always given by Eqs. (l2Ul) . 
(|23p . its behavior at intermediate distances can be qualitatively different, depending on the values of the parameters 
IT and Ifj,. If the temperature is sufficiently low (T < O.lvp/l) and the chemical potential lies in the vicinity of one 
of the Landau levels, the coefficient a in (|22[) is negative and 1/qF *C I <C |o|. In this case the screened potential (TT5)) 
oscillates at intermediate distances 1/qF < r < \a\, as shown in Fig. [3ja),(b). When the chemical potential lies away 
from the Landau levels, or the temperature is larger than OAvp/l, the coefficient a is positive and 4>{r) does not 
oscillate (Fig. [3jc),(d)). In this case the asymptotic behavior of the screened potential for r I is given by 



Ze 
£o 



d 2 q cxp(iqr) 
2tt q + q F + aq 2 2e a(qi - q 2 ) 



[q! [H ( gi r) - Y ( qi r)} - q 2 [H (q 2 r) - Y (q 2 r)} } , (32) 



where q\^ 2 = (1 ± yT — 4aqp)/2a, Hq(z) is the Struve function, and Yq(z) is the Bessel functions of the second kind. 
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Now let us consider the case when both temperature and scattering rate are zero. In this case 

e 2 N f ■ 



11(0,0) 



2ttP 



E ( 2 - 6 °>^ - XM ^ = ^(m) : 



(33) 



n=0 A=± 



where Dq(/i) is the DOS at the Fermi surface for the clean graphene [26|]. When the Fermi level lies between Landau 
levels (which corresponds to the integer fillings) the above expression vanishes, i.e., qp — 0. Restricting ourselves to 
these integer fillings and setting Vt = in (|13[) or, equivalently, setting T = 0, T n = in (fTtj)) . we obtain 



n(o,q) 



E 



)-+ 

' nn 



e 2 N f 



2-Kl 2 ^ M n +M n > 
n.n —I) 



-,(y,A) e 2 N fa , i2 



^ 2 -a 2 )E E E 



:2 „ VM , ^ ^ ^ Mn -\'M n . 

n/=Q n=0 / X'=± 



(34) 



Now the transitions between levels n <-> ±n ± 1 give the main contribution at long wavelengths, because of the 
following asymptotics of the functions (j9|) at y — > : 



0^ +1 (y,A) = Q^ 1 .„( 2/ ,A)=y 



2n + 1 + AA' 



, / nM n+ i 



l)M n 



M, 



n+1 



Q^,(y,A) = 0(y 4 ), n^n±l, Xn ± A' 



This leads to the behavior 



n(0,q)~gaq 2 



iqi « yi ■ 



0(y 4 ), n>0, (35) 

(36) 

(37) 



The coefficient a at zero temperature and scattering rate is always positive and depends on the number of filled 
Landau levels Np and the gap A. It is evaluated to be 



a(N F ,A) 



e 2 N f l 



F(d)+8(p 2 -A 2 )J2 

n=0 



(2-3 0n )(3rc + 2rf) 
\Jn + d 



(38) 
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where d = I A /2v F is the dimensionlcss gap parameter, and we define the function F(d) as 

F{d) = Y(V^+d-Vn-l + df(l + — d =) . (39) 

' \ V^+dVn -l + dj 

At zero gap and Np = 0, we obtain, in agreement with [H, Il8l. W\. 

a( 0) 0) = -^-F(0), ^( )^-6C(-l/2)--^ + O(n- 3 / 2 )^f.247, (40) 



where £(z) is the Riemann zeta function. 

From Eq. (fl8)) we obtain that at long distances the screening is absent (the correction to the bare Coulomb potential 
is of the smaller order), 

<^(r)~— fl-^V r»l. (41) 



V. OTHER LIMITING CASES 



At zero momentum, only the terms with An = X'n' survive in ([7]), and the general expression for the polarization 
function simplifies to 



n=0 \=± 



At zero scattering rate and finite f2 the above expression vanishes, while its limit O — > at finite T n is given by the 
equation ([24| . Therefore, the static long- wavelength polarization function 11(0,0) does not depend on the order of 
taking limits ft — > and q — > 0, unlike in the absence of magnetic field. 

The strong magnetic field limit (I — > 0) of the polarization function ([7]) also depends on the ratio between the 
scattering rate and the frequency. For T n /Q, ^= the main contribution comes only from the lowest Landau level 
(n = 0) and is given by the expression 

\ e 2 N f T f./l \n + \'A + n + iT \ ,/l Xfi + A' A + iT \ 1 

n(n. i) - -2^ n(n + ar ) J^ ± { U + J 2^r J ^(2 + J } ' (43) 

However, this contribution vanishes in the clean graphene limit (more exactly, for To = and nonzero f2). In this 
case the transitions n O — n ± 1 in Q dominate at high magnetic field, resulting in 

II(r!,q)^a(0,0)q 2 , (44) 

Z7T 

which is equivalent to the static long wavelength limit of the polarization function for the clean gapless graphene at 
zero temperature in the case when only the lowest Landau level is filled. 



VI. SUMMARY 



In this paper we have derived the exact analytical expression for the one-loop dynamical polarization function in 
graphene, as a function of wavevector and frequency, at finite chemical potential, temperature, band gap, and taking 
into account the finite scattering rate of Dirac quasiparticles due to the presence of impurities. The most general 
result is given in terms of the digamma function and generalized Laguerrc polynomials and has the form of double 
sum over Landau levels, Eq.Q. In the clean graphene at zero temperature, for the integer fillings of Landau levels, 
this function correctly reproduces the previously obtained results. The derived expression for dynamical polarization 
can be used to calculate the dispersion relation and the decay rate of magnctoplasmons depending on temperature 
and impurity rate. 

The long-range behavior of the screened static Coulomb potential in graphene in magnetic field is found to be 
essentially affected by the presence of impurities or the finite temperature. When either the scattering rate or the 
temperature is nonzero, the usual Thomas- Fermi screening is present, and the resulting potential decays as ~ 1/V 3 , 
which is typical for two-dimensional systems. The strength of the screening oscillates as a function of chemical 
potential or a magnetic field. If both scattering rate and temperature are zero, these oscillations turn into the 
sequence of delta-functions, and for the integer fillings the screening is absent. 
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Appendix A: Calculation of polarization function 

After evaluation of the trace, the equation ([5]) can be written in the following form 



E £ £ 



8w 2 l 4 ~'_ n .Xl J , (iu m + n + iT n sgnu m - \M n )(iu m - s + n + ir n >sgnuj m - s - A'M n >) ' 

71 , Tl — U A, A — zt 77T. — — OO 

(Al) 

where y = q 2 Z 2 /2, and 

= / ^e-^(^) a fflq ,(-^)^(^)^(^) , a = 0,1. (A2) 

The above expression is nonzero only for n, n' > 0. Integrating over the angle and making the change of the variable 
r 2 = 2l 2 t, we get 



I% n ,(y) = 2tt/ 2 / dte-H a ,h(2^t)L a n {t)L a n ,{t) 
Jo 

poo 

= 27rl 2 (-n'-l) a / dte- t J (2 v ^t)L^(t)i-^ Q (t), a = 0, 1, (A3) 



o 



where we have used 



Lt(x) = (-x)- k ^^Lj- + k k (x), l>0, k + l>0. (A4) 



Now, using the formula 7.422.2 in |2i 



OO 

2 i,» / b 2 \ f h 2 

-m+n I " \ Jjs—o 



dx^e-o* J v {bx)L»;°{ax')L° n {ax") = {-l) m+n {2 a y»-^ L°~ m+n [ — )L^ +m ~ n [ — ) , (A5) 



o 

we obtain from (IA3I) 



KAv) = 2rt 2 (-l) n - n (n' + lYe-vL n n - n {y)LlT? a {y) 



2tt1 2 



^±p^e-yy^L^\y)L^:2{y) , a = 0, 1, (A6) 



n>! 



where we again used the formula (|A4[) and the symmetry I" n r(y) = I%> n {y) which follows from (|A2[) . Now we can 
rewrite (|A1[) as 

n(ifi„q) = — ^ 2 E Qnn'(v,m, (A7) 

n.n'— A, A' — ± 

where the functions Q^/(y, A) are defined in (J9j) and 

OO ^ 

j _ \ A (A8) 



m— — oo 
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To evaluate this sum, we expand it in terms of partial fractions and split into four sums in the following way: 



1 



E 



i 



i 



XM n - \'M n > - ifl s - i(T n - r„/) f-f \iuj m + \i + iT n - XM n iu> m -s + M + £T„/ - X'M n > 



1 



XM n - \>M n , - iQ s - <(T n . - r„) 

1 

am„ - A'M n , - i« s - i(r n + r„/) 



m— s 

-1 

E 

771— — OO 

OO OO ■ 

EE 

m— m— s 



i 



l 



iw m + £t — £T n - AM„ iuim-s + M - ilV - A'M„ 
1 1 



(A9) 



iuj m + fi, + iT n - XM n iu> m _ s iT n > - X'M n 



Now, making the change m — > m + s in the first and the last sums, and m — > — m — 1 in the second one, and using 
the summation formula 



( r, _|_ 



1 



n + a n + b 



(A10) 



we obtain 



-TX 



z-, x '>- x (in s ,r,-»,T) 



XM n - x'M n , - in s - i(r n - r„o xM n - x'M n > - ifi s - i(r n , - r„) 
z£(in s , r, t) + z™, (-^^> s , -r, M , r) 

XM n — X'M n ' — ifl s — i(T n + T n >) 
where the functions Z££, (Q, T, fj,, T) are defined in ©. Using the above equation and the relation 



(-<n„ -r, ^ t) = z~? n r* (m s ,r, t) , 

tp(l — Z) = 1p(z) + 7TCOt(7T2:) , 

z^(m s ,r,fi,T) 



-A', -A' , 



which follows from the formula 



we can rewrite (IA7D as 



n(^,q) = ^ E E 0™ : (^ A 



n,n'=0 X.X' = ± 



XM n - X'M n ' - - i(T n - r„ 



z-£'- A (in, J r,- At ,T) 



^(in„ r, /i, t) + z- x 'r y (*n., r, - M , r) 



(All) 

(A12) 
(A13) 

(A14) 



AM„ - A'Af„/ - itt s - i(T n > - T n ) XM n - \'M n > - in s - i{T n + T n ,) 

Making the analytic continuation from Matsubara frequencies by replacing ifl s — > 51 + iO, we finally arrive at ([7]). At 
constant scattering rate r„ = T the result simplifies to 



e 2 N 



n(^)^ E E QtttoX) 



n,n'=0 A,A'=± 



z% (o, r, m, t) + (n, r, r) 

AM„ - A'M„/ - n 

z£(n, r, M , r) + z-*r x ' (n, r, r) 

AM„ - A'M„< - - 2«r 



(A15) 



One can check that the first term in square brackets does not have poles at fl = XM n — X'M n ' since the numerator 
vanishes at this point, 



z£,(n, r, (j,, t) + z-^in, r, t) 



4ir 2 T 



,(\ n-X'M n ,+iT\ ,/l fi-XMn-iT 



2inT 



+ tp' 



2inT 



fi = XM n - X'M nl + e, e -> 0. (A16) 



At r — > the denominators in Eq. (|A15|) become equal, and the overall numerator reads 

z£,(n, o, (i, t) + z~ x ''- x (n, o, r) - ^(n, o, M , r) - z~ x y x ' (n, o, - M , r) 



i 



1 _ fl - XMn 

2 2i7rT 



_ 1 M -AM„ 
1 2 2i7rT 



1 _ /i - X'M n 

2 2i7rT 



_ 1 M~VM n , 
1 2 2inT 



(A17) 



= n F (X'M n >) - n F (XM n ) , 
where we used the property (|A13p of the digamma function. 



11 



Appendix B: Schwinger proper-time calculation of polarization function in magnetic field 



The general expression ([7]) for the polarization function as a double sum over the Landau levels is useful for high 
magnetic hclds. Clearly, for weak fields Eq.([7]) is not convenient since we need to keep many terms in the double sum. 
In general, when L depends on the Landau index n it is impossible even to get a closed expression for the quasiparticlc 
propagator, not to mention the polarization function itself. In principle, it is possible to perform summation in Eq.([3|) 

for L = const and /i 7^ but the expression obtained looks rather cumbersome for further work with it. Therefore, 

00 

we consider in this section only the case T = fj, = 0. Using the identity 1/a = J dte~ at ,a > for introducing the 



proper-time coordinate t, and the formula [3l| 



^LZ{x)z» = (l-z)- a - 1 e X p^ 1 , \z\<l, 



n=0 



we get a closed expression for the fermion propagator: 



S(iu> m ,r) 



dt exp 



f(co 2 



A 2 



4/ 2 



■ coth t 



x I (■yoiUm + A) [P_ (1 + coth t)-P+(l- coth t)] - 

The integrals can be evaluated through confluent hypergeometric functions, 
00 00 
h{a,b) = J dte- at - bcotht = ±e- b r^)y(^,0,2b), I 2 {a,b)=Jdte 



— at — b coth t 



coth t 



dl\ (a, 6) 
db : 



(Bl) 



(B2) 



h{a,b) = 

Hence, we have 

S{iui m ;v 



dt e- at - b coth * coth 2 i - d h (a ' h) ~ - ' ^ 



db 2 



b = 



42 2 ' 



7£ n 

I 2 



-^r(i + |) 



(7 iw m + A) 
* |1- 



2f 



p_r - * -, 1, -5 + p+r 1 + - * 1 + -, 1, 



21' 



a r 

2' '2Z 2 



(B3) 



(B4) 



Using the integral representation (|B2[) for the propagator, we get from © taking the trace and performing the Gauss 
integration over coordinates, 



n(*n a ,cO = — 



e 2 Tl 2 N 



f 



E 



dt dx 



coth t + coth x 



exp 



q 2 l 



2]2 



coth t + coth x 



(A 2 — u; to cj^j)(1 + coth t coth x) + 



Vp (coth t + coth x — q 2 r ) 



I 2 sinh 2 (i + x) 

Introducing new variables, t = z(l + v)/2, x = z(l — v)/2, we obtain 

00 1 



(B5) 



n(*n»,q) 



Til 2 



du 



dv 



exp —it A 



cosh z — cosh zv 



2,2 



2 sinh ; 



(B6) 



sinh z 



where u = l 2 zjv 2 F and the sum 



cosh z — cosh ; 
2 sinh z 



ft 2 

-q^P ) + ucothz[ A 2 + -y 



9 u 9 
9m u dv 



R(u,v,Cl s ) , 



1 — s + sv 



(B7) 
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can be written through the Jacobi elliptic function. For that we use the formula 



m— — oo 



(m+c) 



ic In g 



^ = e lffC "^ 3 (cr|r) = Hr)- 1 ^^! - l/r), q = e 47rT , Imr > 0, (B8) 



where 



6 3 (v, q) = 8 3 (v\t) = 1 + 2 ^ <f 2 cos(27rnw), 



(B9) 



n=l 



and for the third equality we used the Jacobi imaginary transformation. Hence the sum (|B7[) takes the form 



R(u,v,n m ) 



2Ty/¥U, 



Since 



u iL- v iL ]R{ > u > v ^ s) 



1 _ (1 - V)fl m _ 1/(4uT 2 ) 

2 4ttT ' 6 



(1 + i>)» m i/ (4uT 2j 
4ttT 



e -u(l-« a )n;/4 / ! (l + w 2 ) u £} 2 g g 
^ £ Hi 4- u u 

2TJ¥u \ 2 4 5u 



2T,/¥ul 



(BIO) 



(1 -HQOro i/( 4 „ T 2) 

4ttT ' 



(Bll) 



we write 



n(iO s ,q) 



sinh 2 





f dv 


-u ^A 2 


/ 7^ J 


1 y exp 









1 — u cosh 2 — cosh zt> 



1 



cosh z — cosh ; 
2 sinh z 



ucothz A 



2;2 



2 sinh 2 

(l-v 2 )fL 2 s 1 d v d 
4 2u du u dv 



x fl 4 



(1 + l>)£l m i/(4„ T 2 ) 

4T ' 



(B12) 



The above integral is divergent at u = reflecting the primitive divergence of the polarization function. Therefore, in 
order to get finite result one should regularize the initial expression, for example, by subtracting the same expression 
with A replaced by M — ¥ oo (the Pauli-Villars rcgularization) which means that we write 



II(zf2 s , q) = lim 



e 2 Nf f du f dv 



m^oo 2tt 3 / 2 1 2 j v^j 2 



(A 2 M 2 ) 



(B13) 



Carefully separating the part with M 2 and taking into account that 



lim 



exp(-uM z ) — + M- 
1 2u 



1 

2u 



= 0, 



(B14) 



we finally get the following expression for the polarization function at finite temperature in a magnetic field, 

OO 1 s 

e 2 N f f du f dv J exp(-uA 2 ) 



n(ifi.,q) 



2ir 3 / 2 l 2 J 2 

o -1 ^ 



sinh , 



zexp 



1 — v 2 „ cosh z — cosh zv 
tt 2 + — — q 2 v 2 F 



sinh; 



cosh 2 — cosh zv 



q 2 l 2 I + coshz 



2 sinh ; 



/2A 2 Z 2 Q 2 J 2 



\ Vp 2v 2 F sinh2z 



2z sinh z 
coshz coshzu — 1 



2 sinh 2 



x0 4 



-1/(4«T 2 ) 

4nT ' 



cosh 26*4 



0ie -l/(4«T 2 ) 



(B15) 



where we also performed the integration in parts of terms with derivatives over u, v. 
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Now we consider several limiting cases of Eq. (|B15[) and compare them with expressions existing in the literature. 
Taking the limit T — > is very easy since theta-functions turn into units. After some transformations the zero 
temperature limit can be recast in the form 



n(zft s ,q) 



e 2 Nfcf f du f dv zcosh 
4tt 3 / 2 



' coth z sinh ; 



x exp 



A 2 + 



sinhz 

cosh z — cosh zv 
2z sinh z 



q 2 4 



the result first obtained in Ref . [27j . 

On the other hand, taking the limit of zero field, I — > oo, in Eq. (|B15[ ) we get 



n(in„q) 



e 2 N f 



2tt 3 / 2 v% J u 3 / 2 



du f dv 
2 



exp 



2 + u 2A 



X 4 



(l+v)tts l/(4«T 2 ) 



4ttT 



, e 



0, 



0)e -l/(4«T=) 



e-« A - 1 



(B16) 



Q 2 S + v 2 q 2 v 2 F 



(B17) 



The integration over u in (|B17[) can be performed explicitly using a series representation for thcta functions, we get 
in terms of the integration variable x = (1 + i>)/2: 



IT^q) 



e 2 N f 
2ttv f 



dx 



n 2 + q 2 v 2 F + 4[A 2 - x{l - x)q 2 v 2 F ] sinh(a(x)/T) + ^ cosh(A/2T) 



4a(x) 



£>(!) 



2L>(a;) 



(B18) 



where 



a(x) = ^A 2 + x(l - x)(Q 2 + q 2 v 2 F ), D(x) = cosh 2 (a(x)/2T) - sin 2 (7rsx). 



This expression can be rewritten in somewhat different form if we integrate by parts the last term in square brackets 
and then use the identity among the integrals, 



ATQ S J dx\a[AD{x)} = 2VL S J dx 



a(x) smh(a(x)/T) 



+ (n 2 s +q 2 v 2 F ) I dx(l-2x) 



sin(27rsa;) 
D(x) ' 





which can be obtained following the method described in the appendix A of Ref. [32[ ■ Finally, we have 

l 



n(»n a ,q) 



2tt n 2 + q 2 v 2 F 



dx 



2Tlog[4D(a;)] 



A 2 sinh(o(x)/T) 



a(x) D(x) 



(B19) 



(B20) 



For A = 0, Eq. (|B20| is in agreement with Eq.(A20) (together with (A23), (A26)) in jH] while for T = it reduces 
to the well known expression for the vacuum polarization operator in QED3 [33j. 
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